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ABSTRACT 

We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21 . 

These images were obtained with the Subaru HiCIAO camera, adaptive optics, and the polarized differential imaging 
technique. We resolve the disk in scattered light at H - band for stellocentric OH < r ^ Of! 6 (12<r<75 AU). We 
compare our results with previously published spatially resolved 880 pcm continuum Submillimeter Array images 
that show an inner r < 36 AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor 
invoked to explain SR21 ’s sub-mm cavity cannot be “universal” for all grain sizes. Even significantly more moderate 
depletions ( 8 = 0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (5 ~ 10 -6 ) 
are inconsistent with our //-band images when they are assumed to carry over to small grains, suggesting that surface 
grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. 

In fact, the radial polarized intensity profile of our H - band observations is smooth and steeply inwardly-increasing 
(r -3 ), with no evidence of a break at the 36 AU sub-mm cavity wall. We hypothesize that this profile is dominated 
by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the 
previously postulated existence of a sub-stellar companion to SR21 at r ~ 10-20 AU, and find that we can neither 
exclude nor verify this scenario. This study demonstrates the power of multi wavelength imaging of transitional 
disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for 
clearing these disks of their large midplane grains. 

Key words: planetary systems - protoplanetary disks - stars: individual (SR21) - stars: pre-main sequence - 
techniques: polarimetric 
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1. INTRODUCTION 

Transitional disks are of particular interest to the astronomi- 
cal community because they have features intermediate between 
optically thick gas-rich protoplanetary disks and optically thin 
dusty debris disks, suggesting that they may be an intermediate 
evolutionary step. Transitional disks were first classified based 
on their spectral energy distributions (SEDs) as objects with 
significant mid-infrared excess emission indicative of starlight 
reprocessed by a circumstellar disk, but a lack of strong near 
infrared (NIR, ~l-5 /xm) excess. This NIR deficit implies a 
dearth of the very warmest grains in the inner disk (Strom 
et al. 1989; Skrutskie et al. 1990). spatially resolved sub-mm 
imaging of a number of transitional disks has revealed that 
many of them are indeed cleared of large thermally-emitting 
dust grains out to 10AU or more, well beyond the dust sub- 
limation radius, thereby suggesting that an additional clearing 
mechanism may be at work (Andrews et al. 2009; Hughes et al. 
2009; Brown et al. 2009; Andrews et al. 2011, referred to as 
All throughout the remainder of this paper). 

The precise mechanism by which these clearings are made is 
widely debated in the literature, with several viable alternatives 
put forward, including clearing by massive forming protoplanets 
(e.g., Dodson-Robinson & Salyk 2011; Zhu et al. 2011), 
clearing by stellar mass companions (e.g., Ireland & Kraus 
2008; Jensen & Mathieu 1997), photoevaporation (e.g., Owen 
et al. 2012; Clarke et al. 2001; Hollenbach et al. 1994), and 
grain growth (e.g., Birnstiel et al. 2012; Dullemond & Dominik 
2005). Each theoretical explanation for the observed presence 
of sub-mm disk cavities makes its own predictions for clearing 
timescales, efficiencies and radii. However, the small amount of 
available spatially resolved multiwavelength data on transitional 
disks has limited the ability of modelers to test these theories 
against observation. As such data accumulate, they reveal 
an increasingly complex picture of the “cleared” regions. In 
particular, cavities seen at long wavelengths do not necessarily 
manifest themselves in NIR scattered light (e.g., this work, 
Dong et al. 2012). Spectroscopic studies also reveal a variety of 
gaseous species in the inner regions of transitional disks (Salyk 
et al. 2009; Pontoppidan et al. 2008). As large dust grains, small 
dust grains, and gas all appear to respond quite differently to 
the mechanism or mechanisms that create transitional disks, a 
thorough understanding of their clearings can only be grasped 
through accumulation of a statistically significant sample of 
multiwavelength data on transitional disks. 

Thus far, the most success in spatially resolving transitional 
disks has been achieved in the sub-mm regime (e.g., Al 1), with 
a small number of additional studies using mid-infrared (MIR, 
10-40 /xm) aperture masking (e.g., Eisner et al. 2009, hereafter 
E09), NIR spectroastrometry (e.g., Pontoppidan et al. 2008, 
hereafter P08) and MIR interferometry (e.g., Ratzka et al. 2007) 
to infer inner disk structure. Efforts are currently underway to 
image a sample of transitional disks at shorter wavelengths, 
including the Strategic Exploration of Exoplanets and Disks 
with Subaru (SEEDS; Tamura 2009) campaign. 

Despite initial success at resolving the rim of the LkCa 15 
cavity in scattered light (Thalmann et al. 2010) with a Locally 
Optimized Combination of Images (LOCI) method, additional 
transitional disks imaged by SEEDS are revealing that these 
cavities, including that of LkCa 15, are more elusive when 
viewed in polarized light (Hashimoto et al. 2012; Muto et al. 

* Based on data collected at Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan. 


2012; Dong et al. 2012; Mayama et al. 2012; Grady et al. 
2013; J. Wisniewski et al. 2013, in preparation). On the surface, 
this suggests that NIR scattering material is present, and even 
abundant, inside the cavities. Although sub-mm studies are 
generally revealing that transitional objects are indeed cleared or 
depleted of large grain material inside of the cavity radii inferred 
from their SEDs, NIR “clearings” in transitional disks may not 
be as ubiquitous nor as fully cleared as their sub-mm cousins. 
In fact, it appears that some transitional disks contain abundant 
small grain material inside of observed sub-mm cavities. 

Polarized NIR scattered light and sub-mm imaging, though 
both sensitive to circumstellar dust, probe very different regimes 
in the disk, both in terms of density and height above the mid- 
plane. The sub-mm emission from transitional disks is ther- 
mal, tracing mainly mm-sized grains that are concentrated 
(“settled”) toward the disk midplane (e.g., Dullemond & 
Dominik 2004; Brauer et al. 2008; D’ Alessio et al. 2006). These 
large (a > 10 /xm) grains are very important for understanding 
the large scale structure of the disk, as they house most of the 
disk’s dust mass and their coagulation drives the process of 
planet formation under the core accretion scenario (e.g., Pol- 
lack et al. 1996). The density and temperature of the grains 
at the disk midplane also determine the size and shape of the 
remainder of the disk. 

At NIR wavelengths, however, it is unclear whether scattered 
light emission originates from an optically thick small grain 
disk surface, or from more diffuse material in or above the disk. 
This uncertainty that is compounded by the small number of 
spatially resolved disks in the literature. It is unknown whether 
disks with large grain clearings in the sub-mm are also cleared 
of small grains, nor is it understood what proportion of the 
disk mass lies in the small grain population. The fact that we see 
axisymmetric polarization structure in the NIR from transitional 
objects reveals that they house circumstellar disks with some 
small grain material (Whitney & Hartmann 1992). However, 
very little dust mass is needed to scatter at H- band, so it is 
unclear how much mass resides in this portion of the disk and 
whether it is optically thin or optically thick in this wavelength 
regime. If transitional disks are optically thick in the NIR, as is 
the case for protoplanetary disks, then this scattered emission 
will probe only the t = 1 surface. Small grains are believed 
to be vertically extended by a factor of two or more relative to 
larger grains (Dullemond & Dominik 2004; Brauer et al. 2008), 
so NIR scattering in an optically thick disk should occur well 
above the disk midplane. It is unclear whether or not clearing 
mechanisms acting on the large grains at the disk midplane 
would have an observable effect on the NIR emission under this 
scenario. 

Conversely, if transitional disks are optically thin in this 
wavelength regime, as is the case for debris disks, then scattered 
light emission probes the entire vertical extent of the disk. In 
this case, scattered light emission originates from small grains 
throughout the vertical extent of the disk, down to and including 
the disk midplane. Observations of the same disk at both 
sub-mm and NIR wavelengths are compelling in that they allow 
for comparison of the settled and upper surface layers of the disk, 
the large (mm-sized) and small (/xm-sized) grains, separately. 
A good disk model should thus be able to self-consistently 
reproduce observations at both wavelengths simultaneously, as 
well as any additional multiwavelength data. Our efforts to do 
just that for the SR21 disk are described in detail in Section 4. 
Section 2 describes our data collection methodology, and our 
result are summarized in Section 3. 
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1.1. Observational History ofSR21 

SR21A is a 2.5 Af 0 , ~1 Myr old Class II young stellar object 
(YSO) in the Ophiuchus Star Forming Region. Recent astromet- 
ric observations suggest that the distance to the Ophiuchus star 
forming region is 1 19 zb 6 pc (Lombardi et al. 2008). However, 
we adopt a distance of 1 35 pc (Mamajek 2008) in order to remain 
consistent with previous studies. SR21 has been classified as a 
6" 7 binary (Barsony et al. 2003); however, Prato et al. (2003) 
concluded that the two components were not coeval, and that the 
B component does not exhibit infrared excess emission and is 
therefore unlikely to host a circumstellar disk of its own. The A 
and B components also have discrepant proper motions (Roeser 
et al. 2010), so the “binary” does not appear to be a bound, co- 
eval system. We will therefore refer to SR21A as simply SR21 
throughout the remainder of this paper. 

SR21 was first identified as a transitional disk by Brown et al. 
(2007), who inferred a 0.45-18 AU gap in the disk based on a 
fit to the ~ 1-100 /I m SED. In order to explain the NIR excess 
seen at <5 /xm, they invoked a small uncleared region in the 
inner disk from 0.25 to 0.45 AU. Since this initial identification, 
a variety of puzzling and seemingly contradictory data have 
emerged that continue to favor a disk that is not simply cleared 
of all material in this 0.45-18 AU region. 

The first evidence of a lack of complete clearing in the inner 
disk came from CO rovibrational spectroastrometry done with 
CRIRES on the Very Large Telescope by P08, who observed 
molecular gas emission originating from a region well inside 
of the dust gap inferred by Brown et al. They concluded that 
this emission arose from a narrow 7-7.5 AU ring of molecular 
gas. It is not surprising to have found gas emission inside of 
the dust gap, as SR21 may still be accreting, albeit weakly 
(cKT 8 - 84 Mq yr -1 , a relatively high upper limit; Natta et al. 
2006). However, the physical mechanism causing an apparently 
truncated ring of molecular gas at this radius is still an open 
question. P08 suggest that the inner 7 AU cutoff corresponds 
to a physical truncation of the disk, while a rapid drop in 
gas temperature outward of 7 AU, perhaps due to a thick 
disk component in the region, causes the CO flux to quickly 
drop below the detection threshold. This explains the observed 
truncation of CO emission at 7.5 AU without requiring that the 
disk material be physically truncated at this radius. 

E09 also postulated that there is emission from within the 
~18AU “cavity” based on MIR (8.8 and 11.6 /xm) aperture 
masking data taken with TReCS on Gemini South. They put 
forward an alternate interpretation to Brown et al. (2007), 
capable of fitting both the SR21 SED and their visibility 
data. They suggest that the disk is completely cleared within 
10 AU, but with a warm companion at 10-18 AU whose 
circum(sub)stellar material creates the NIR excess seen in 
the SED. 

The first spatially resolved imaging study of SR21 was done 
at 340 GHz (880 /xm) and was published by Brown et al. (2009) 
among a sample of four disks with large central clearings. 
Follow-up data, also at 340 GHz, were presented in Andrews 
et al. (2009). We base our sub-mm modeling efforts on All, a 
detailed imaging and modeling study of 12 transitional disks, 
including SR21. Using two-dimensional Monte Carlo radiative 
transfer modeling to fit the sub-mm data and the infrared SED, 
All concluded that the 36 AU sub-mm central cavity in SR21 
was depleted in material by a factor of MO -6 relative to the 
outer disk. 

As multiwavelength data accumulate on SR21, it is clear that 
its structure, particularly in the region interior to the 36AU 


sub-mm clearing, cannot be understood through SED modeling 
alone because there are several degenerate solutions capable of 
explaining it (Brown et al. 2007; Eisner et al. 2009). Spatially 
resolved images are needed at multiple wavelengths, probing 
multiple regions of the disk, before good theoretical constraints 
can be put on the mechanism clearing the interior of SR2 1 . In the 
next section, we describe new spatially resolved NIR scattered 
light imaging, which contributes further to the puzzle of the 
SR21 “cavity.” 

2. OBSERVATIONS AND DATA REDUCTION 

Polarized differential imaging (PDI) of SR21 was done in 
the H band (1.6 /xm) on 2011 May 22 with the high-contrast 
imager HiCIAO (Tamura et al. 2006; Hodapp et al. 2006) and 
adaptive optics system AO 188 (Minowa et al. 2010) on the 
8.2 m Subaru telescope on Mauna Kea. These observations 
were conducted as part of the larger SEEDS survey, which was 
begun in 2009 October and is currently entering its third year. 

The data were taken in “qPDI” mode in which a double 
Wollaston prism is used to split the beam into four 512 x 
512 pixel channels (0 / . / 0095 pixel -1 ), two each of o- and 
e-polarizations. The splitting of each polarization state into 
two separate channels reduces saturation effects on the Hawaii 
2RG detector. In order to obtain full polarization coverage and 
minimize artifacts, a half waveplate was rotated to four different 
angular positions (0°, 45°, 22?5 and 67?5). This cycle was 
repeated 18 times on SR21, with a 15 s exposure per waveplate 
position, for a total of 18 minutes of integration time. 

Each image was bias subtracted, flat fielded and bad pixels 
removed in the standard manner for infrared data analysis using 
custom IDL scripts. Pinhole data were then used to create a 
distortion solution for each channel, which was applied before 
the images were rotated to a common on-sky geometry. All of 
the 15 s science exposures of SR21 were saturated to a radius 
of 6-8 pixels. For this reason, all subsequent combinations have 
had an inner r — 8 pixel mask applied. 

Each individual channel was then aligned using Fourier cross- 
correlation. These aligned images were combined to create 
Stokes Q and U images using standard differential polarimetry 
methods, namely adding together each set of two identically 
polarized channels and then subtracting these combinations 
from one another. The four channels were also added directly 
together to form a total intensity image. 

To estimate the performance of the Subaru AO system during 
our observations, we used two unsaturated images of the point- 
spread function (PSF) reference star HD 148212, and 12 and 
10 unsaturated images of SR21 taken before and after the main 
SR21 data set respectively. The average Strehl ratio of the HD 
148212 images is 0.27, and the average Strehl ratios of the 
SR21 before and after images are 0.22 and 0.20 respectively. 
The standard deviations of the Strehl ratios in the two SR21 
data sets are 0.018 and 0.016 respectively, indicating that AO 
correction was maintained steadily. 

2.1. Polarized Intensity Images and Polarization Vectors 

The polarized intensity (PI) image shown in Figure 1 is a 
median combination of individual PI images for each of the 1 8 
waveplate cycles, computed via PI = y/ Q 2 + U 2 . Polarization 
angles were calculated for each image according to the formula 
Op = 0.5 tan -1 (U / Q) and median combined. These are shown 
overplotted on the PI image in five pixel bins in the lefthand 
panel of Figure 2. The length of each vector has been scaled 
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Figure 1. Median polarized intensity (PI) image of SR21 at H- band. The gap radius inferred by All and the companion orbit inferred by E09 are indicated with 
dashed ellipses. These are circles inclined at 14° and rotated to a major axis PA of 100° (All) and 60° (E09) (indicated with the solid white and red lines). North is 
up; east is left. The central (y/07 were saturated in the raw images and have been masked out. A diffraction limited spot is shown at the bottom right for reference. 
Modeling suggests that the emission deficit ~0?2 south of the central source is a result of the disk geometry and is not a physical structure. 

(A color version of this figure is available in the online journal.) 
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Figure 2. Left: polarization vectors overplotted on the median polarized intensity (PI) map of SR21. A general centrosymmetric geometry is evident in the SR21 
vectors, suggesting that this emission does indeed trace scattered light from a circumstellar disk. Right: polarization vectors for the PSF star HD 148212 overplotted 
on its median PI image, which has been scaled by a factor of three to compensate for the difference in exposure time. In this case, the vectors do not show a coherent 
pattern, and PI emission is weak and compact. This is interpreted as photon noise residuals and is used to estimate the error in the polarized intensity profiles shown 
in Figure 3. In both cases, vectors are scaled by the strength of the polarized intensity signal in the region and are not plotted when they fall at or below noise level. 
Vectors whose lengths correspond to 500, 1000, 1500 and 2000 mJy arcsec -2 are shown in the bottom right of the SR21 panel as a scale reference. 

(A color version of this figure is available in the online journal.) 


according to the magnitude of the PI in that bin. A first 
order instrumental polarization correction has been applied, as 
described in detail in Hashimoto et al. (201 1). 

As has been demonstrated by other groups (Close et al. 1997 ; 
Oppenheimer et al. 2008; Apai et al. 2004; Quanz et al. 2011), 
PI images are a powerful data product for studying circumstellar 
disks. Unlike most other varieties of disk observations, PI 


data are generally spared the difficult and imperfect task of 
PSF subtraction under the assumption that direct starlight is 
randomly polarized, while light scattered off of a circumstellar 
disk shows a preferred polarization direction according to 
the scattering location in the disk. Specifically, the scattering 
process will preferentially scatter light polarized perpendicular 
to the direction to the illuminating source, while light coming 
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Figure 3. SR21 polarized intensity image (background) and polarization vectors (overplotted) after first order correction for a seeing limited halo component. 
This method reveals the polarization vectors to be more centrosymmetric, particularly along the disk minor axis, than originally inferred. Local deviations from 
centrosymmetry persist in the northeast and southwest quadrants of the disk, where unresolved disk structures and/or companions could induce local polarization 
structure that deviates from the larger centrosymmetric pattern. The process of subtracting this halo component has introduced additional PSF artifacts into the 
background PI map, most evident in the negative astigmatic pattern centered on the star. We have thus chosen to model the uncorrected polarized intensity map, which 
we believe more closely resembles the true structure of the SR21 disk. Vectors whose lengths correspond to 500, 1000, 1500 and 2000 mJy arcsec -2 are shown in the 
bottom right of the SR21 panel as a scale reference. 

(A color version of this figure is available in the online journal.) 


directly from the star will be unpolarized and randomly oriented. 
Extended polarized emission is therefore unlikely to occur in 
systems without scattered light disks, and light scattered from 
a disk should show a centrosymmetric pattern of polarization 
vectors. 

Although direct starlight is randomly polarized when it 
reaches the telescope, multiple reflective surfaces in the HiCIAO 
optical path are capable of polarizing some fraction of this 
incoming light. This phenomenon is referred to as instrumental 
polarization (IP). In order to estimate the quality of IP removal 
via standard techniques in our data, we reduced our data on the 
PSF star HD 148212 using the same methodology as SR21, 
including a first-order instrumental polarization correction. 
HD 148212 is a G2V star at d = 105 ± 15 pc and -10' 
from SR21. It was observed just before SR21 on the same 
night and in the same observing mode, with 17 cycles of 
waveplate rotations and 5 s per exposure, for a total of 5.7 
minutes of exposure time. As shown in the righthand panel of 
Figure 2, no polarized emission was observed beyond r > O'.' 2 
in the initial first order IP corrected HD 148212 PI image. 
Although the polarization signal inside of this radius is small 
compared to the PI of SR21 in the same region, it does show the 
characteristic four-lobed astigmatic pattern of the HiCIAO PSF. 
Additionally, while polarization vectors for HD 148212 do not 
show centrosymmetry, they do show some alignment along the 
vertical axis of the instrument. These are both hints that there is 
a residual uncorrected instrumental polarization signature, and 
further correction is warranted. 

To first order, we can approximate the residual instrumental 
polarization by measuring the average polarization strength and 
direction of the signal in the HD 148212 PI image. We assume 
that HD 148212 is an ordinary diskless unpolarized star, and 
that all of its extended emission is due to residual instrumental 
polarization or photon noise. We construct an artificial “halo” 
by multiplying the average strength (1.1%) in the HD 148212 
halo into the raw intensity image, then creating artificial “halo 
Q” and “halo U” images from it using the measured average 
polarization angle (—11°). We combine these synthetic Q and 


U images to create a “halo PI” image and subtract it from the 
original HD 148212 PI image. We also subtract the “halo Q” 
and “halo U” images from the raw Q and U images, and use 
the resulting images to calculate corrected polarization vectors, 
which we overplot on the halo- subtracted PI image for HD 
148212 in the righthand panel of Figure 2. The extended PI 
emission from the PSF star is reduced to r < O'.' 1 through this 
process, and the resulting vectors are randomized, no longer 
showing a preferred polarization direction. We are therefore 
confident that this provides a reasonable approximation of the 
residual instrumental polarization. 

The need to apply this same “halo removal” procedure to the 
SR21 PI data can be seen in the Figure 2 SR21 polarization map. 
Although the SR21 vector pattern is clearly centrosymmetric 
along and around the disk major axis even in the uncorrected 
image, it deviates somewhat from centrosymmetry along the 
same instrumental axis as the aligned vectors in the HD 148212 
image (approximately the minor axis of the disk). Since not all 
SR21 polarized emission is spurious, approximating the strength 
and direction of the halo is nontrivial in this case. We used a 
region with low intrinsic disk emission to construct a “halo” for 
subtraction, namely 100 pixel (0795) wide regions centered on 
the minor axis at 30 pixels < r < 50 pixels ((X'285 < r < 07475) 
on either side of the disk. Following the same procedure as with 
HD 148212, we created PI “halo” Q, U and PI images from the 
polarization direction (—8? 8) and strength (1.02%) measured in 
this region and subtracted. 

Subtraction of this first-order approximation reveals a more 
centrosymmetric pattern, as evidenced in the lefthand panel of 
Figure 3. However, it also introduces significant PSF artifacts 
to the PI map, as can be seen in the background PI contours of 
Figure 3. For this reason, we have chosen to model the “raw” 
unsubtracted radial profile in Section 4, noting, however, that 
some of the emission interior to O'.' 15 is likely to be contaminated 
by a small degree of instrumentally polarized stellar flux. The 
polarized halo strength of 1.02% is an order of magnitude 
smaller than our lower limit for the intrinsic SR21 polarization 
(10%), so this contamination is minimal. 
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SR21 Intensity 
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Figure 4. Intensity images of SR21 (left) and HD 148212 (right). HD 148212 has been scaled according to differences in exposure time and intrinsic H - band magnitude 
(based on published 2MASS values) to match SR21 for subtraction. Some uncorrected astigmatism is evident in the four lobes of the PSF in both stars. 

(A color version of this figure is available in the online journal.) 


Despite its crudeness, the “halo subtraction” method high- 
lights an interesting feature of the SR21 polarization vector 
structure, which is that several regions in the disk with high 
signal (notably the northeast and southwest quadrants) show 
deviations from centrosymmetry even after halo subtraction. 
This may indicate that there is unresolved structure in the disk 
at these locations. Interestingly, the deviation in the northeast is 
consistent with the location of the companion inferred by E09, 
which could be explained if the companion had an accretion disk 
of its own, with its own polarization structure superimposed on 
that of the SR21 disk. 

2.2. Intensity Images 

Intensity images were also constructed from the PDI observa- 
tions of SR21 by cross correlating, centroiding and aligning the 
images in each of the four channels at each of the four waveplate 
(0°, 45°, 22?5 and 67?5) positions for each of the 18 observa- 
tion cycles. The four channels at each waveplate position were 
summed for each cycle, and the resulting 1 8 images were then 
median combined to create an SR21 intensity image. 

Unlike PI images, total intensity images are dominated 
by stellar flux, and PSF subtraction is necessary. For this 
purpose, we created intensity images of HD 148212 using 
the same methodology described above. The final reduced 
(but unsubtracted) median combinations of the SR21 and 
HD 148212 total intensity images are shown in Figure 4. Both 
reveal significant uncorrected astigmatism, as evidenced by 
the four-lobed pattern centered on the star. The SR21 disk is 
not visible in the unsubtracted total intensity image. The PSF 
subtraction process is described in detail in Section 3.3. 

3. RESUFTS AND ANAFYSIS 

Our principle observational results are the final polarized 
intensity (Figure 1) and total intensity (Figure 4) images and 
polarization vectors (Figures 2 and 3). These images were 
analyzed in order to estimate the extent, inclination, major axis 
orientation and polarization percentage of the SR21 disk. The 
derivation of each of these properties is described in detail in 
this section. 


3.1 . Polarized Intensity Isophotes 

We began by fitting elliptical isophotes to our PI data in order 
to provide our own measure of the inclination and major axis of 
the SR21 disk, independent of the values in the literature, which 
were determined at other wavelengths. These elliptical fits are 
shown overplotted on Figure 5 in white, beginning at a radius of 
10 pixels from the center (well outside the saturation radius) and 
extending outward to 50 pixels, where the PI flux approaches the 
noise level. Errors in isophotal fits were determined by summing 
in quadrature the intrinsic error estimation in the IRAF 29 ellipse 
procedure, which was used to fit the isophotes, and the standard 
deviation of the parameter across the image. 

The average ellipticity of these isophotes suggests a disk 
inclination of 14 zb 2° (where 0° would be a face-on geometry). 
This is smaller than the value reported by All; however, their 
estimate of 22° was made with the caveat that the sub-mm 
ring seen in SR21 is very narrow, and the geometry is not 
well constrained by their observations. The H - band emission is 
significantly more spatially extended; therefore, we adopt 14° as 
the preferred value in our modeling efforts, which are described 
in Section 4. There is no statistically significant deviation in our 
calculated inclination value between the inner and outer disk 
of SR21. 

There is a clear counterclockwise progression from the inner 
to the outer disk in the major axis of the elliptical isophotal fits 
(overplotted in white in Figure 5). The average orientation of the 
major axis over the entire disk is 86° zb 1 1 ° . The large uncertainty 
is due to the fact that the PA migrates from inner to outer 
disk. This happens relatively smoothly, but with a pronounced 
jump between 17 and 20 pixels. Inside of this radius, which is 
of particular interest because it corresponds to the location of 
E09’s hypothesized companion, the average major axis position 
angle is 76° zb 8° (measured east of north), while outside of 
this radius it is 96° zb 7°. This shift in the average major axis 
PA of ~20° between these two regions, and of 42° in total 
between inner (68°) and outermost (110°) isophotes, suggests 
that the SR21 disk may contain unresolved structures such as 
a warp or a spiral arm. Some of this PA migration might be 

29 IRAF is distributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science Foundation. 
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SR21 Polarized Intensity Isophotes 
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Figure 5. Elliptical isophotal fits to the SR21 polarized intensity image. The fits are shown overplotted on the PI image in white starting at a semimajor axis of 10 pixels 
(well outside the saturation radius) to a semimajor axis of 50 pixels (where the PI flux drops to noise level). The major axis of each isophote is overplotted as a straight 
white line. The average value of the major axis position angle is clocked by ~15° from the inner to outer disk, suggesting unresolved disk structure. The ellipticity of 
the isophotes was also used to estimate the inclination value of 14° (where 0° would be face-on). 

(A color version of this figure is available in the online journal.) 


Table 1 

Position Angle Estimates 


PA 

n 

Radius Probed 
(AU) 

Reference 

15 ±4 

7.0 ± 0.4 

P08 (v = 1-0) 

16 ±3 

7.6 ± 0.4 

P08 (v = 2-1) 

44 ±7 

11 ± 1 

E09 

61 ± 8 

15 ± 1 

E09 

76 ±8 

12-25 

This work 

96 ±7 

25-65 

This work 

100 

30-120 

All 

105 

30-120 

Brown et al. (2009) 


explained by instrumentally polarized emission with the shape 
of the HiCIAO PSF, particularly because the most pronounced 
turnover occurs at approximately the same radius as the extent 
of the HD 148212 halo. However, the observed SR21 isophote 
major axis migration continues well beyond that radius, and we 
therefore conclude that it is a real feature of the data. 

We combine our results with previous estimates of the major 
axis PA from the literature in Table 1. The All estimate for 
the major axis PA of 100° (overplotted on Figure 5 in red) is 
consistent with our measurement of the PA of the outer disk. We 
find that previous estimates of the PA of the disk are strongly 
dependent on the radius being probed, and our observation of 
a counterclockwise trend with radius fits well into the broader 
picture. This suggests that there is unresolved non-axisymmetric 
structure in the SR21 disk, and provides further confirmation 
that the trend we observe is not caused purely by residual 
PSF-shaped instrumental polarization signal. 

3.2. Radial Brightness Profiles and Asymmetries 

We used the final median combination of unsaturated total 
intensity images and the published Two Micron All Sky Survey 


(2MASS) photometry values of SR21 and HD 148212 to es- 
tablish a flux conversion for our data. In doing so, we make 
the assumption that the H - band magnitudes of both objects re- 
mained steady in the decade between 2MASS and HiCIAO 
observations. H - band emission of T Tauri stars has been 
shown to be variable (Carpenter et al. 2002); however, the 
ADU-to- Jansky conversion factor calculated using the 
HD 148212 images differs by only 1.6% from that calculated 
using the SR21 images; therefore, we believe that this is a robust 
assumption and have adopted the average of the two values as 
our flux conversion factor. We have applied a software mask to 
the saturated region interior to r = 8 pixels; however, astigmatic 
PSF artifacts extend out to r ~ 13 pixels (O'/ 13, ~18 AU). Our 
data are therefore most robust at 18 AU < r < 80 AU. 

Computed flux-converted radial PI profiles for SR21 are 
shown in Figure 6. Brightness profiles are shown for both the 
azimuthally averaged case (filled black points) and for the case 
where we limit analysis to a region 11 pixels wide centered 
on the disk major axis (open red points). We also overplot 
(open blue points) the radial profile of a PSF-subtracted PI 
image. Although PSF subtraction is nominally unnecessary for 
PI images, our analysis in Section 2.1 revealed that residual 
instrumental polarization contributes at a low level. We therefore 
utilized the PI image of the PSF star HD 148212, and scaled it 
such that the noise in the background region was equivalent to 
that of SR21 and subtracted. This provided a subtraction with 
fewer evident PSF artifacts than our earlier “halo subtracted” 
PI image and serves to verify that the slope and steepness of 
the radial PI profile is insensitive to contamination by residual 
instrumental polarization. 

In all three cases, we plotted the profile in five pixel radial bins 
excluding the innermost two bins (so interior to O'/ 1), which have 
pixels inside the saturation radius. This binning was chosen to 
match the diffraction limit of (X/05 at the H band. The brightness 
profile along the major axis is systematically higher than the 
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SR 21 Radial Polarized Intensity Profiles 



radius in arcseconds 


Figure 6. Radial polarized intensity profiles for SR21 binned to the H - band diffraction limit (075) and shown as triangles. The black filled triangles represent the 
azimuthally averaged radial profile and the red open triangles represent the profile as measured in an 1 1 pixel wide region centered on the major axis of the disk. The 
blue open triangles represent the profile of the PSF-subtracted image. All three profiles follow a very similar slope, with the major axis profile systematically higher in 
flux per pixel by 20% relative to the azimuthally averaged profile, and the PSF-subtracted image systematically lower, as is to be expected. Because the images were 
saturated interior to 0707, radial profiles are only shown for bins that fall completely outside of the saturated region (from 07 1 outward). Also shown is a canonical r~ 2 
profile in green and ar -3 profile in blue. The r -3 profile is a better fit to the data, although there is some evidence that it becomes shallower in the outer disk beyond 
0745. The veracity of this shallower profile in the outer disk is brought into question by the steepness of the PSF-subtracted profile, suggesting that the flattening in the 
outer profile may be due to halo effects. Errors were estimated from the photon noise in the four individual channels used to create the PI image summed in quadrature 
(black error bars). The location of the sub-mm cavity wall at 36 AU and the noise floor are indicated with vertical and horizontal dashed lines, respectively. 

(A color version of this figure is available in the online journal.) 


azimuthally averaged profile, as is to be expected based on 
the increased brightness of the disk in this region; however, 
the shapes are virtually identical. The PSF-subtracted profile is 
systematically lower than that of the unsubtracted profile, as is to 
be expected. Unlike the other two profiles, the PSF-subtracted 
profile does not flatten in the outer disk, suggesting that this 
trend may be caused by residual instrumental polarization in 
the unsubtracted profiles. We have elected to fit the azimuthally 
averaged profile of the unsubtracted image in our modeling, 
which is described in Section 4. 

These radial PI profiles reveal that H- band scattered light 
emission exists in the SR21 disk out to a distance of at least 
~076, which is comparable to the disk extent in the All 
Submillimeter Array (SMA) images. This corresponds to a 
physical scale of ~80AU, or more than twice the radius of 
the observed sub-mm cavity. The excess in the NE quadrant 
of the disk, which could be due to the contribution of a 
companion as discussed in Section 3.3, is outside of the region 
used to create the major axis radial profile, and is smoothed 
out in the azimuthally averaged profile, therefore does not 
have any noticeable effect on either of our derived radial 
profiles. 

Errors in the radial PI measurements were estimated using 
both the background photon noise in the four individual qPDI 
channels of the SR21 raw images and using the observations of 
the PSF star HD 148212. The residual PI flux surrounding HD 
148212 was summed in the same manner as the SR21 radial 
profile and is comparable to the photon noise. The error bars in 


Figure 6 and all of the radial profile plots shown in the modeling 
section correspond to the photon noise. 

A canonical r~ 2 power law for the surface brightness of a 
flared disk is shown overplotted in green on all radial profiles 
and is a poor fit to the data. The best fit radial profile to the 
data is of the form r -3 . Although there is some indication 
that it may flatten beyond 0745, this trend is not exhibited in 
the PSF-subtracted PI profile, so we believe it to be spurious. 
The smooth, steep incline in PI from the outer disk in to the 
saturation radius is the crux of the new NIR data, and explaining 
it is the main driver of the modeling efforts described in 
Section 4. 

3.3. Intensity Images and Compatibility with the 
Existence of a Substellar Companion 

Direct addition of the four qPDI channels yields a total 
intensity image, which may also be diagnostic of the disk, but 
is dominated by stellar emission. PSF subtraction is required 
in order to isolate the disk component in these images. PSF 
subtraction is an imperfect process in systems such as SR21 for 
a number of reasons. First, the Strehl ratio may vary between 
the two sets of observations, and within each set itself, and 
Strehl ratio mismatch can introduce artificial structure into the 
PSF-subtraction. The SR21 and HD 148212 data sets were 
calculated to have comparable Strehl ratios (~20%); therefore, 
we do not believe that Strehl ratio mismatch is a dominant effect. 
Secondly, speckles that survive the median combination process 
are not necessarily cospatial in the two data sets; therefore, 
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SR21 PSF Subtracted Intensity 
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Figure 7. H - band PSF-subtracted intensity image with the 18 AU orbit of the companion proposed by E09 overplotted. There is a distinct excess of emission in the 
northeast quadrant (A) at a location that is consistent with their hypothesized companion, although its proximity to PSF artifacts makes robust identification impossible 
with this data. There is also a northward extension of the disk that is nominally present in the PI image (B). The image has been smoothed with a FWHM = 6 pixel 
Gaussian in order to suppress artificial structures induced by speckle noise. 

(A color version of this figure is available in the online journal.) 


local excesses and deficits must be viewed with some suspicion. 
Thirdly, the SR21 data set was centrally saturated, while the 
HD 148212 data set was not, and so the counts in the SR21 image 
are an underestimate of the total emission. Finally, and perhaps 
most importantly, the total intensity of SR21 includes both a 
direct stellar and scattered light disk contribution, while HD 
1482 12’s emission corresponds to direct starlight alone. These 
last two effects make estimation of the appropriate scaling of the 
PSF image for subtraction somewhat ad hoc, but by employing 
a series of reasonable approximations we can place bounds on 
the appropriate scale factor for subtraction. 

We chose to scale the PSF star HD 148212 first by the 
relative ratio between its 2MASS H - band magnitude and that of 
SR21, a quantity which we will call a. Assuming that the total 
H - band magnitude of these systems has not varied significantly 
with time, subtracting the PSF star intensity image scaled by 
this value is an over subtraction, both because some portion of 
SR21’s H - band magnitude corresponds to its disk contribution 
and because the SR21 image is saturated at the center. Therefore 
we use this subtraction (SR21 — a*HD 148212) as an “upper 
limit” on the disk’s total intensity in scattered light. 

We then created a grid of PSF subtractions with successively 
lower scale factors applied to the PSF star. We used the ratio 
of scattered polarized light to total scattered light in the disk as 
a metric in establishing the lower limit for the PSF scale fac- 
tor. This ratio of polarized to total scattered light ( P = PI/ I& S k) 
should be on the order of 50% in the case of 100% efficient scat- 
tering of unpolarized starlight off of spherical grains; however, 
variables such as porosity, oblateness, and alignment as well 
as particle size and composition effects, multiple scattering ef- 
fects, wavelength dependence, and inclination effects have been 
invoked to explain the observed deviations from this value (e.g., 
Whitney & Wolff 2002; Simpson et al. 2009). Observed scat- 
tered polarization fractions vary from ~10% to ~70% in tran- 
sitional disks, with typical values of approximately 40% (Dean 


Hines and Glenn Schneider, private communication based on 
HST program GO 10178, 2012 October). 

We deem values of <10% to be unphysically low based on 
these results; therefore, we apply the scale factor resulting in 
P = 10% (0.5a) as our lower limit for PSF subtraction. Shown 
in Figure 7 is the PSF subtraction resulting from a scale factor 
0.75a, which is halfway between the upper (a) and lower (0.5a) 
limits established above. The image has been smoothed by a 
Gaussian with FWHM = 6 pixels to mitigate the effects of 
speckles on the inferred structure of the disk. 

In estimating the ratio of polarized to total disk intensity 
P (PI/ /disk), we use this intermediate PSF scaling and an annulus 
of 15 < r < 60 pixels. This region begins just outside the radius 
of the PSF astigmatism and ends at approximately the outer 
radius of the PI disk. Within this region, we have rejected all 
pixels where (1) the PI flux is below noise level, (2) P > 1.0, 
indicating a PI speckle, or (3) P < 0, indicating an I speckle. 
Our upper and lower limits for the subtraction tell us that SR21’s 
P value lies somewhere between 10% and 27%, and a = 0.75 
gives P = 18%. 

The final PSF-subtracted and masked intensity images of 
SR21 for all scale factors investigated reveal two structures, both 
of which are also present in the PI data. The first is an excess 
in the NE quadrant near the location of the E09 companion, 
and the second is an arc extending upward toward the north. 
These are marked as A and B in Figure 7. Although it is entirely 
possible that these features are spurious, they are intriguing in 
the context of the E09 observations and our own observation of 
the trend in PA with radius in the SR21 disk, which is suggestive 
of structure. 

E09 inferred a warm (~700 K), extended (~40 R Q ) com- 
panion 18AU from SR21 based on the system’s mid-IR flux 
and visibilities. Such a source would contribute ~5% of SR21’s 
bolometric flux, or ~0.4 L©, and would need to be a low-mass 
star surrounded by a dusty shroud or disk. The enshrouded 
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companion would likely be considerably brighter in the H band 
than a ~700 K blackbody, itself a magnitude 16 source at SR21’s 
distance. We note that a bolometric luminosity of 0.4 L© is char- 
acteristic of a ~0.2 Mq T-Tauri star at 1 Myr, which would have 
an effective temperature of ~3000 K (Kenyon & Hartmann 
1995). The observed excess at point A in our data could also be 
due to a simple overdensity in this region of the disk. 

We obtained our images of SR21 with the image rotator off, 
allowing us to use Angular Differential Imaging (ADI) to search 
for a companion. A traditional ADI analysis (Marois et al. 2006) 
detected no statistically significant sources. Because of SR21’s 
declination, we were only able to obtain ~10° of field rotation, 
making it difficult to apply more sensitive techniques like 
Locally Optimized Combination of Images (LOCI; Lafreniere 
et al. 2007). We therefore attempted a variation of LOCI using 
the ACORNS-ADI pipeline (Brandt et al. 2013), in which 
we supplemented the comparison frames with 30 principal 
component images culled from a large library of observations 
of single stars. 

We detect no statistically significant companion candidates 
down to a limiting contrast of 100 at O'/ 14 (18 AU at the 
distance to SR21). This is not a meaningful limit for a substellar 
companion. However, it does likely rule out a stellar companion 
sufficiently luminous to account for the mid-IR excess reported 
by E09. Follow-up observations at longer wavelengths or with 
extreme adaptive optics will be necessary to reach fainter 
detection limits. 

4. MODELING 

To assist in interpreting our results, we ran a suite of disk 
models using the three dimensional radiative transfer Whitney 
code (see Whitney et al. 2003a, 2003b, for details). The key 
capabilities, modifications from previously reported versions, 
and input assumptions of the Whitney code are described in 
Section 4.1 

In Section 4.2, we describe our efforts to reproduce the 
sub-mm observations of All using an independent code, and 
our verification of the Al 1 best fit parameters for the large grain 
dust disk. 

In Section 4.3, we report on several degenerate explanations 
for the three features that make the SED of SR21 unique: (1) 
the steep 12-20 /rm rise, (2) the marked NIR excess, and (3) the 
very small size of the 10 fim silicate feature. 

In Section 4.4 we report on our attempts to reproduce the 
observed steep r -3 polarized intensity radial profile, and its 
ability to constrain the geometry and optical depth of the disk. 

Section 4.4 also details our efforts to merge these three pieces 
into a single self-consistent model, and in Section 5 (Discussion) 
we speculate on the implications of our best fit model. 

4.1. The Whitney Code and Input Assumptions 

The Whitney code mimics settling and grain growth by speci- 
fying the composition and structure of a “large grain” disk sepa- 
rately from that of a “small grain” disk. The large grain disk can 
be artificially settled relative to the small grain disk by setting 
its scale height parameters to lower values. Although both disks 
contain a distribution of grain sizes, we will refer to them simply 
as the “large grain disk” and “small grain disk” throughout. 

The Whitney code treats both thermal emission and scattered 
light in the disk. Grain temperatures are computed using 
the method of Lucy (1999) and the code treats anisotropic 
scattering using a Henyey-Greenstein function (Henyey & 


Greenstein 1941). The albedo, average scattering angle, and 
maximum polarization are pre-computed at each wavelength 
based on a Mie calculation for the ensemble of grains in each 
population. The polarization phase function is approximately 
dipolar, following White (1979). 

The Whitney code utilized herein has been modified from 
previously reported versions and is described in detail in Dong 
et al. (2012) and Whitney et al. (2013). The modifications 
of particular relevance to this paper were threefold. First, the 
disk density distribution was modified in order to allow for 
an exponential cutoff in mass as a function of radius, making 
it identical to the disk mass distribution of models run by 
All. Secondly, the new models allow for a “puffing up” 
of the inner disk edge and the inner edge of the disk gap, 
a physical situation widely invoked to explain NIR spectral 
excess (e.g., Dullemond & Monnier 2010) and increasingly 
invoked to explain the observed features of gapped disks (e.g., 
Jang-Condell & Turner 2012). The final modification made to 
the code allows for the differential depletion of large and small 
grains, in addition to independent geometries between the two 
disks (a feature included in earlier versions of the code). This 
means that the disk radial density exponent (a, p a r~ a ), radial 
scale height exponent (/3, Z a r - ^), absolute scale height 
at 100 AU, and depletion (5) are all specified independently 
inside and outside of a disk gap or cavity. In effect, this allows 
for four independent disk components — large and small grain 
disks inside and outside of a specified gap radius. The size 
distribution and composition (and therefore opacity as a function 
of wavelength) of large and small grains are also specified 
independently. 

The basic model input parameters (e.g., stellar mass, radius, 
temperature) were determined based on values in the literature, 
and those adopted by All in particular. The choices of model 
parameters are given in Table 2. In total, over 250 independent 
models of the SR21 disk were calculated in this modeling effort. 

Initial input assumptions included a disk mass of 0.006 M 0 , 
a gas: dust ratio of 100:1, and a large grain dust mass fraction 
of 85%. Neither of these values is well constrained by our 
H - band observations, so we have assumed the same values as 
All. We have fixed the composition of the settled dust disk to 
the grain prescription of Wood et al. (2002), a power-law grain 
size distribution with an exponential cutoff (Kim et al. 1994) 
composed of amorphous carbon and astronomical silicates with 
solar abundance constraints for carbon and silicon (C/H~3.5 x 
10 4 and Si/H ~ 3.6 x 10 5 , respectively; Anders & Grevesse 
1989; Grevesse & Noels 1993). The functional form of this size 
distribution is 



where a is the grain size, n is the number of grains, and the 
parameters p, a c , and q, which control the distribution shape, are 
adjusted to fit the wavelength dependence of the observations. 
In particular, the long wavelength slope of the SED and the peak 
of the MIR excess are sensitive to these parameters. Following 
the Wood et al. (2002) best fit model for the HH30 disk, we use 
p = 3.0, q = 0.6 and a c = 50 /xm. This is a good match to the 
grain properties inferred by Al 1. 

Several small grain dust prescriptions were also investigated, 
and the observed SED and H - band images were found to be 
relatively insensitive to their particulars. We use the inter- 
stellar medium (ISM) dust prescription of Kim et al. (1994) 
which has a smaller average grain size than our other models. 
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Table 2 

The Range of Whitney Model Input Parameters Investigated for SR21 


Model Input Parameters 

Parameter 

Model Input(s) 

Input Assumptions 

37s tar 

2 Mq* 

^star 

3.15 Rq* 

Ltar 

5750 K + 

Ay 

6.3 

AJdisk 

0.006 Mq* 

r gap 

None, 18, 36* AU 

r disk, in 

r S ub - 7 AU 

r disk, out 

200 AU 

i 

14, 22* 

rc 

15 AU* 

Large Grain Disk Parameters 

Composition 

Wood et al. (2002)+ 

fd 

0.85* 

Zioo 

1.52 AU* 

^outer 

2.15+ 

Pouter 

1.15+ 

a gap 

2.15+ 

^gap 

1.15+ 

5 

10" 8 , 10 -7 , 10 -6 \ 10" 5 , 10" 4 , 10" 3 , 10 -2 , 10" 1 , 1 

Small Grain Disk Parameters 

Composition 

Kim et al. (1994); Cotera et al. (2001) 

fd 

0.15 

Zioo 

2-15 AU 

^outer 

0.5-4 

Pouter 

1-3.15 

^gap 

-5-10 

ft gap 

1-3.15 

8 

10“ 8 , 10“ 7 , 10“ 6 , 10“ 5 , 10“ 4 , 10“ 3 , 10“ 2 , 10 _1 , 1, 2, 5, 10 


Notes. Where several discreet values were attempted, the table entry is given in the form of a list, and for parameters 
where a range of values was attempted, the table entry gives the full range of attempted values. Parameters identical 
to those of Al 1 are marked with an asterisk, while equivalent parameters (rounded for convenience, compatibility, 
or translated to Whitney code syntax) are marked with a + sign. Rc is the characteristic radius for the exponential 
density cutoff as defined in All ,/j is the fraction of the disk’s dust mass, Zioo is the disk scale height at 100 AU, 
a is the exponent for the radial midplane density dependence, p is the exponent for the radial scale height 
dependence, and <5 is the depletion factor of the gapped region relative to the outer disk. 


Like the large grain disk, the small grains follow a power law 
size distribution with an exponential cutoff; however, in this 
case they are composed of silicates and graphite, with maxi- 
mum grain sizes of 0.16 /zm and 0.28 /zm respectively. 

4.2. Sub-mm Image Modeling 

Our H - band PDI observations probe only scattered light from 
the upper surface of the disk and therefore do not constrain the 
properties of the large grain disk where the majority of the 
disk mass lies. However, we have endeavored to match 
the sub-mm observations of All in addition to our newly- 
acquired data in order to develop the most physically motivated 
disk model possible. As a first test of the compatibility of the Al 1 
models with the Whitney code, we began with input parameters 
for the disk that were identical (or nearly) to the best-fit model 
of All. 

Figure 8 shows the SMA image taken by All in panel 1, 
and the equivalent 880 /zm image produced by our model in 
panel 2. We are able to produce a similarly well-fitting model 
with the Whitney code and identical input parameters. Model 
outputs were sampled with the same spatial frequencies ob- 


served with the SMA, and images were created with the same 
inversion and deconvolution parameters as the observations. 

Having shown that the All and Whitney models produce 
equivalent results with identical input assumptions, we then 
investigated the effect of gap depletion on the sub-mm images to 
test the robustness of the Al 1 assumption that the gap is depleted 
by a factor of 8 = 10 -5,9 relative to the outer disk. Keeping all 
other disk parameters fixed, we varied only the depletion factor 
8 in the interior of the gap by an order of magnitude for each 
model from 8 = 10 -8 (two orders of magnitude higher depletion 
than the All assumption of ~10 -6 ) to zero depletion (which, 
at least on the surface, is the most consistent with the evidently 
undepleted cavity seen in our observations). 

We break these models into three categories: “high” depletion 
models ( 8 = 10 -8 -10 -6 ), “moderate” depletion models ( 8 = 
10 -5 -10 -3 ) and “low” depletion models (8 = 10 -2 -10°) 
and show the real 880 /zm image of All, our simulated 
image, residuals, and visibility fits for each of these models 
in Figures 9-11. In each case, the model images were scaled 
such that the total flux was equivalent to that of the real SMA 
observation, and the required scaling is given in the model panel 
of the figures. 
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Figure 8. Left: the 880 /xm SMA image of All. Right: a large grain disk model with All’s best-fit input assumptions (5 = 10 -6 , M disk 0.006 Af 0 , r gap = 36 AU, 
etc.) produced by the Whitney code. 

(A color version of this figure is available in the online journal.) 



The “low” depletion models provide poor fits to both the 
880 pm image and the visibility profile. We conclude, as did 
All, that the SR21 disk is depleted in large sub-mm emitting 
grains interior to 36 AU. All of the moderate and high depletion 
models provide reasonable fits to both the 880 pm image and 
the visibility profile; however, the 8 = 10 -5 model shows the 
lowest residuals and the 8 = 10 -6 model provides the best fit 
to the visibility profile. The 8 = 10 -6 model requires just a 
2% increase in flux to sum to the observed flux density (easily 
explained by small differences in input opacity assumptions), 
while the 8 = 10 -5 model requires an 8% scaling, so we adopt 
the 8 = 10 -6 value for the large grain r < 36 AU disk in all 
further modeling. This value is closest to the All best fit value 
of 8 = 10 -5,9 ; therefore, we consider our large grain models to 
be functionally equivalent. 

As a test of the robustness of the large grain disk assumptions, 
we verified that the small grain disk parameters specified in 
our models did not create observable variation in the simulated 
880 pm image. As expected, we found that wide variation in 
input small grain parameters has no observable effect on the 
880 pm images produced by our models. The 880 p m emission 
of the large grain disk is therefore strongly decoupled from the 
small grain disk parameters. 

Having verified that we can reproduce the All result, that 
the All inferred gap radius and depletion factor for large 
sub-mm emitting grains is a good match to the models produced 
by both codes, and that this result is unaffected by variation 
in small grain parameters, we treat the All large grain disk 
parameters (with the exception of disk inclination, for which we 
adopt our measurement) given in Table 2 as fixed throughout 
the remainder of our modeling efforts. 

4.3. SED Modeling 

There are sufficient data in the literature to construct a well- 
sampled observed SED for SR21. In addition to 2MASS /, 


H, and K photometry, SR21 photometry has been obtained 
with Spitzef s Infrared Array Camera (IRAC) and Multiband 
Imaging Photometer (MIPS) at 3.6, 4.5, 5.8, 8 and 70 pm, 
the AKARI spacecraft at 18 pm, the Infrared Astronomical 
Satellite (IRAS) at 12, 25 and 60 pm, the SMA at 880 pm, and 
the James Clerk Maxwell Telescope’s Submillimeter Common- 
User Bolometric Array (SCUBA) at 350, 450, 850 and 1300 pm. 
10-35 pm spectra from the Spitzer Infrared Spectrograph (IRS) 
and 52-97 pm from the MIPS SED mode are also available. 
These are shown overplotted on all SEDs generated by our 
models. Although we show it for completeness, the MIPS 70 pm 
photometry point has been excluded from our fits due to its 
inconsistency with both the MIPS SED spectrum and IRAS 
60 pm photometry. All reported photometric points and spectra 
have been extincted according to a standard ry = 3.1 interstellar 
extinction law (Cardelli et al. 1989). The extinction for SR21 
has been estimated to be as high as A v = 9.0 (Prato et al. 
2003), but All have argued that the published photometry is 
most consistent with a lower value of Ay ~ 6.3. We also find 
that lower extinctions are most consistent with the observed 
optical/NIR SED, so we adopt the same A v = 6.3 value 
for SR21. 

The SR21 SED shows circumstellar dust excess emission 
longward of the stellar blackbody peak, as evidenced by the 
long wavelength excess (see Figure 12). Unlike protoplanetary 
systems, which tend to have flatter transitions between the stellar 
and dust peaks in the SED indicating significant dust populations 
at all temperatures, the SR21 SED has a distinct double peaked 
profile, reaching a local minimum near 12 pm. This suggests 
that there is a dearth of hot grains in the disk, a conclusion 
supported by the observed sub-mm cavity. Debris disk systems, 
on the other hand, are much more depleted of small grains, and 
their emission traces the stellar photosphere out to wavelengths 
of 20 pm or more. In contrast, the SR21 SED rises rapidly 
beyond 12 pm, indicating a healthy population of warm grains 
that emit at X > 12 pm. 
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Figure 9. Highly depleted (8 = 10 _8 -10 -6 ) large grain disk models produced by the Whitney Code. From left to right, the panels are (a) the observed SR21 SMA 
image of Al 1 , (b) the final Whitney model image, where the raw model output was sampled with the same spatial frequencies observed with the SMA, (c) the residuals, 
and (d) the model fit (in red) to the SMA visibilities. Input parameters were identical to those of Al 1, except that the depletion factor (5) for large grains inside of 
36 AU was varied. The depletion factor is given above each model row. A scale factor was also applied to each model image to equalize the model and SMA fluxes 
before subtraction from the observed image. This factor is given in panel (b) of each model. The minor discrepancies between modeled and observed fluxes can be 
easily explained by slight differences in input opacities. The bottom image is the same best fit model shown in Figure 7. 

(A color version of this figure is available in the online journal.) 


Shortward of the 12 /xm local minimum in the SR21 SED, 
the emission is also greater than photospheric, which indicates 
a small population of very warm grains (a hot inner disk 
component) or, alternatively, a warm companion. The SR21 
SED shows little evidence of a 10 /xm silicate feature, while 
many similar transitional disk systems exhibit strong, broad 
silicate emission at this wavelength. A good disk model for 
SR21 therefore needs to reproduce all three of these distinct 
SED features: (1) the steep 12-20 /xm rise, (2) the sub- 12 /xm 
excess, and (3) the lack of a strong, broad 10 /xm silicate feature. 

The Whitney code produces stable SEDs with relatively few 
(~ 10 6 ) photons, making it much less computationally expensive 
than generating images/radial profiles. We were therefore able 
to explore much of the available parameter space for SR21 to 
find several classes of reasonable models capable of reproducing 
the observed features of the SED. 


4.3.1. NIR Excess Fits 

The NIR excess shortward of 10 /xm is a distinct feature in the 
observed SED. The (small) excess in this region was originally 
attributed to a 0.07-0.25 AU warm inner disk component. 
Indeed, adding this component to our models produces a 
well-fitting NIR excess. However, the MIR observations of 
Pontoppidan et al. (2008) are somewhat inconsistent with the 
presence of this interior hot disk component. Their data show a 
sharp truncation in gaseous CO at 7 AU, suggesting a high flux 
of UV photons. This feature would be difficult to produce with 
an inner hot disk component at r < 7 AU capable of intercepting 
them. Furthermore, allowing the large grains to fill the cavity 
in to 0.07 AU in our models, even while depleting them by a 
factor of 10 -6 , results in point- like sub-mm emission at the disk 
center, which is not present in the All observations. 
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Figure 10. Moderately depleted (8 = 10 5 -10 3 ) large grain disk models produced by the Whitney Code. See Figure 9 caption for details. 

(A color version of this figure is available in the online journal.) 


The presence of the warm companion predicted by E09 is 
another potential mechanism for creating excess NIR flux in 
SR21, and in this case allows for a larger inner small grain 
disk truncation radius. Our observations can neither confirm 
nor exclude the presence of a companion, yet they do provide 
indirect support of its existence, as detailed in Section 3.3. 

Otherwise identical models with an (1) 7 AU truncation and a 
warm (T = 730 K, R = 30 r Q ; E09) blackbody companion (2) 
a disk that extends inward to the sublimation radius at 0.07 AU 
are shown in Figure 12 to demonstrate that both assumptions 
reproduce the observed NIR excess. Without any means of 
probing disk structures interior to the saturation radius of 0? 1 
in our observations, we cannot break the degeneracy between 
inner disk extent and the existence of a companion in the NIR 
SED fits using scattered light imagery. 

We hold a slight preference for the warm companion hypoth- 
esis based on the P08 CO results and our own observations of 
excess intensity and deviation from polarization axisymmetry 
(indirectly supporting the companion hypothesis and described 
in detail in Section 3.3), so we assume a 7 AU truncation in 


all of our scattered light models. We add the warm companion 
component to all of our modeled SEDs in order to reproduce the 
observed NIR excess, but acknowledge that a significant inner 
disk extent is also a plausible mechanism for creating this fea- 
ture. We do not add PI emission from the companion into our 
simulated images because the PI contribution of such a com- 
panion is (1) not well constrained by models and (2) should 
not have a significant effect on the azimuthally averaged radial 
profile. 

4.3.2. 10 iim Silicate Fits 

The final feature of note in the SR21 SED is the lack of a 
prominent 10 /xm amorphous silicate feature, which tends to 
be rather strong and broad in Herbig disks and is believed to 
originate from the surface layers of the inner (<20AU) disk 
(van Boekel et al. 2005; Meeus et al. 2001). The IRS spectrum 
spans this region of the SED and also that of the weaker 18 /xm 
silicate feature without any evidence of a strong line flux at 
either wavelength. 
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Figure 11. Minimally depleted ( 8 = 10 2 -10°) large grain disk models produced by the Whitney Code. See Figure 9 caption for details. 

(A color version of this figure is available in the online journal.) 


We find that the majority of the 10 /xm silicate emission comes 
from the inner disk wall in our models, and we can suppress it 
by either restricting the scale height of the small grain inner disk 
wall or by restricting small grains to a very small proportion of 
the total mass (< 1%). Generally speaking, many otherwise well- 
fitting models overproduce 10 /xm silicate emission. Changing 
the small grain prescription may be warranted in future modeling 
attempts. Specifically, the silicate features could be suppressed 
by either increasing the grain sizes or limiting the proportion 
of silicate grains in the disk. We leave this for future modeling 
studies. 

Our scattered light imagery cannot give any information about 
the size of the small grain disk wall where the 10 /xm feature 
originates, whether it lies at 7 AU or 0.07 AU (the sublimation 
radius). We limit the size of the inner disk wall to small scale 
heights in all of our scattered light models in order to suppress 
the 10 /xm feature, but in some cases not severely enough, 
resulting in overproduction of 10 /xm emission in the SEDs. 
We could presumably correct this by running additional models 
with smaller inner wall scale heights or reduced small grain 


populations, but have elected not to focus on reproducing this 
feature precisely for models that are poor fits to the radial PI 
profile due to the computational expense. The quality of fit to 
the 10 /xm feature was taken into account in choosing the final 
best- fitting disk model, described in Section 4.4 and shown in 
Figures 15 and 16. 

4.3.3. Fits to 10-20 fim Rise 

The steepness of the rise from 12 to 25/xm is perhaps the 
most unique and significant feature of the SR21 SED and is the 
primary feature that led to its “transitional disk” classification. 
The presence of this feature requires that there be a sharp 
increase in disk emitting area (e.g., a “wall”) at intermediate 
(20-40 AU) radii. However, the precise geometry of the “wall,” 
as well as its composition, is not well constrained by the 
observed SED. 

The simplest way to create a “wall” is to allow for depletion 
inside of the r — 36 AU gap, so that the disk makes a transition 
from optically thin to optically thick at this radius. We know 
that the disk is optically thin at r < 36 AU for long wavelengths 
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Figure 12. Model NIR SEDs produced by otherwise identical models with and without (bottom and top panels respectively) the addition of a 700 K blackbody 
companion to the SED and with either a 7 AU inner disk truncation (right) or a disk extending inward to the sublimation radius (left). These reveal that the k <12 /xm 
NIR excess can be reproduced either by extending the disk inward to the sublimation radius (upper left) or by truncating the disk at a larger radius and adding in the 
NIR contribution of a 700 K blackbody companion (lower right). 

(A color version of this figure is available in the online journal.) 


(and therefore in large, thermally-emitting grains), and optically 
thick beyond; therefore, the large grains certainly do present a 
“wall” at 36 AU. All of the models that we consider are therefore 
“gapped” in large grains. 

However, the sub-mm data places no such constraint on small 
grain depletion. In fact, there is reason to believe that small grain 
dust filtration occurs at gap radii (Zhu et al. 2012; Rice et al. 
2006), and that material must survive into observed sub-mm 
dust gaps, as a number of gapped disks are still observed to 
be actively accreting (Najita et al. 2007; Espaillat et al. 2012; 
Salyk et al. 2009). The observed H - band radial profile of the 
SR21 disk would seem to support this hypothesis, as it shows 
no evidence of a change in small grain disk properties at r = 
36 AU. Until our Subaru H - band data were obtained, there were 
no spatially resolved data at shorter wavelengths to constrain 
the properties of the small grains and the amount to which they 
are depleted inside of the sub-mm cavity. 

The simplest hypothesis, put forward by previous SED 
modelers (Brown et al. 2007, All, E09) is that the small 
grain disk follows the same structure as the large grain disk, 
including the r < 36 AU gap. This hypothesis reproduces 
the observed SED, but it is inconsistent with our new Subaru 
H - band PI profile. Simulated radial profiles and SEDs for 
various small grain depletion factors are shown in Figure 13. 


They show that although small grain depletion in the r < 36 AU 
region reproduces the 12-20 fim rise in the SED, it is strongly 
inconsistent with the observed H - band radial profile. All models 
in which the inner disk is strongly depleted, and therefore 
optically thin, suggest that we should have seen an excess in 
H - band emission at the cavity wall and a turnover in the radial 
profile at this radius. 

Undepleted and low depletion models, on the other hand, are 
inconsistent with the observed SED. Additionally, their radial 
profiles are too shallow to be consistent with the H - band data. 
The Appendix details our efforts to find a scenario in which 
the small grain disk, whether optically thin or optically thick at 
r < 36 AU, contributes to the 12-20 /xm rise in the SED. 

In the end we find that the steep 12-20 /xm rise can be created 
with a disk depleted only in large grains in its interior. Our data 
suggest that the small grain disk is not “gapped.” 

4.4. Scattered Light Imagery 

At first glance, the H - band radial PI profile is difficult to 
reconcile with the large, heavily depleted sub-mm cavity, as 
it spans the region of the sub-mm wall without any evident 
break or turnover. In fact, the principal features of the observed 
radial PI profile, namely its smoothness and its steepness 
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Figure 13. Top: model SEDs of varying r < 36 AU depletion factors overplotted on the observed SR21 SED. Five “Gapped” model runs of different depletion factors 
(<5) are shown with solid colored lines, ranging from 10 -8 depletion in red to zero depletion in purple. In the case of the heavily depleted <5 = 10 -8 and <5 = 10 -6 
models, the NIR excess in the SED can be reproduced either by extending the disk inward to the sublimation radius or by adding the contribution of a 700 K blackbody 
companion. This is discussed in detail in Section 4.3.1 Bottom: H - band radial polarized intensity profiles for the same models overplotted on the observed profile. The 
flux of the undepleted model was fixed to the data at r = 36 AU, and all other models scaled by the same amount, in order to highlight changes in the radial profile at 
the cavity wall. None of the models provide a good fit to the data, revealing that the small grain disk is not only undepleted, its geometry is different from that of the 
large grain disk. 

(A color version of this figure is available in the online journal.) 


(~r -3 ), were not reproducible under any of the scenarios we 
investigated in which the small grain disk contained a depleted 
cavity. Furthermore, we were unable to reproduce the radial 
profile under any scenario in which the H - band scattered light 
originates from the surface of an optically thick small grain 
component, whether gapped or not. This effort is described 
in detail in the Appendix, and representative models for each 
variety of optically thick disk model investigated are shown 
in Figures 17-19. This is a surprising result, as we had no 
reason to expect that transitional disks are optically thin at short 
wavelengths. 

We find that the steepness of the radial profile requires that 
we invoke an optically thin small grain component that spans 


the entire disk. In order to create this, we must either (1) relax 
the All assumption that the small grains compose 15% of the 
disk’s dust mass and restrict their contribution to <1%, making 
the small grain dust disk optically thin throughout, or (2) invoke 
a third optically thin disk atmosphere or envelope component 
that has a larger vertical extent than either the large or small 
grain disks. 

Under scenario (1) the optically thick “wall” producing the 
12-20 fim SED rise must be composed entirely of large grains 
in order to allow the small grain disk to be optically thin. 
As we noted in Section 4.3, we do not require a small grain 
wall in order to reproduce the 12-20 /xm rise in the SED. 
In fact, we can reproduce the SED structure at A, > 12 /xm 


17 


The Astrophysical Journal, 767:10 (26pp), 2013 April 10 


Follette et al. 


Small Grain Optically Thin Disk Model 


Large grain density distribution 


Spectral Energy Distribution 



Small grain density distribution Radial Polarized Intensity Profile 



radius in arcseconds 


Figure 14. Sample model outputs for a representative optically thin disk model. The lefthand panels represent the density of large grains (upper) and small grains 
(lower). The large grain model is identical to the All best-fit model described in Section 4.2 except that the dust wall at 36 AU has been inflated by a factor of two 
in order to account for the lack of small grain contribution to the 12-20 /zm rise in the SED. The small grain model has an inner truncation at 7 AU. The righthand 
panels provide model fits to the observed SED (upper) and H - band radial polarized intensity profile (lower). In this case, the lack of 10 /zm silicate emission is due 
to the very small amount of small grain material populating the inner disk. As revealed by the radial profile, an optically thick small grain disk alone cannot mask the 
discontinuity at the 36 AU large grain dust wall. As justified in the text and shown in Figure 15, an underlying optically thick small grain layer is needed to mask this 
discontinuity, while this optically thin layer is required to reproduce the r~ 3 radial polarized intensity profile. 

(A color version of this figure is available in the online journal.) 


with a large grain gap alone. We do need to increase the 
scale height of the gap wall in order to provide enough flux 
at these wavelengths (relative to models where an optically 
thick small grain component contributes to the rise), but not by 
such a degree as to lead us to declare these models unphysical 
(Z W all> 1AU). 

Since an optically thin small grain disk means that we in- 
evitably “see” the large grains at the midplane, these mod- 
els require that the large grain disk not scatter efficiently at 
the H band. Otherwise, they suffer from the same local max- 
imum problem that led us to exclude “Gapped” small grain 
disks. As can be seen in Figure 14, which shows a scenario 
under which the small grain disk is optically thin throughout, 
the Wood et al. (2002) dust prescription chosen for the large 
grain disk (and found to be a good match to that of All) does 
scatter weakly at the H band, resulting again in a turnover in the 
radial profile at 36 AU. These models would require an alternate 
grain size distribution and/or composition in order to prevent a 
local excess at 36 AU. We have not endeavored to explore this 
possibility here. 


We disfavor explanation (1) because it requires that small 
grains be optically thin throughout the vertical extent of the 
disk. An entirely optically thin small grain disk is very unlikely 
given the small amount of disk mass allowed and the observed 
optically thick r > 36 AU large grain component. In order to 
be optically thin throughout, we require not only that <0.1% of 
the disk mass lie in the small grain disk, representing a likely 
unphysical degree of grain growth or small grain depletion, but 
also that the large grain disk’s grain size distribution contain no 
grains that scatter at the H band, an equally unlikely scenario. 

Under scenario (2), which we believe to be far more likely, 
an envelope component lies above the optically thick small 
grain dust disk with a significant vertical extent. This envelope 
contains very little total mass (on the order of 10 -7 M 0 ) 
and is composed of small ISM-like grains (Kim et al. 1994) 
Representative model outputs for this scenario are shown in 
Figure 15, and a schematic of the disk structure is shown in 
Figure 16. 

Although we call this third disk component an “envelope,” 
it is not clear that it fits into the canonical definition. 
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Figure 15. Sample model outputs for a representative three component disk model. The lefthand panels represent the density of large grains (upper) and small grains 
(lower). The large grain model is identical to the All best-fit model described in Section 4.2. The small grain disk consists of two components: (1) an underlying 
optically thick small grain disk of the “Steeply Curved” variety with a 7 AU inner truncation and (2) an overlying vertically-extended optically thin component that 
contributes thirty times more scattered light than the optically thick component. The righthand panels provide model fits to the observed SED (upper) and H - band 
radial polarized intensity profile (lower). The presence of the optically thick small grain disk acts to mask the discontinuity at the large grain dust wall that was revealed 
in Figure 14. In this case, the 10 /zm silicate feature is small due to the small vertical extent of the optically thick small grain disk wall and the small amount of material 
populating the more vertically extended optically thin small grain component. 

(A color version of this figure is available in the online journal.) 


There has been no evidence in the literature to suggest that 
SR21 is still actively accreting, and the only constraints on the 
accretion rate are upper limits (Natta et al. 2006). This is in con- 
trast to a number of other gapped disks that do show evidence 
of accretion (Salyk et al. 2009). 

It seems that this optically thin component may be more 
readily explained as a vertically extended disk “atmosphere” of 
sorts, rather than an accretion envelope. Whether this is part 
of the natural structure of the disk, a remnant of the envelope, 
or a disk wind launched from the optically thick disk surface 
cannot be determined from our data. In fact, our observations are 
relatively insensitive to its precise geometry. The radial profile 
of the optically thin component follows an r -3 power law in 
all cases where the disk a — /3 = — 1 , which allows for a large 
range in flaring and density structures. 

Under this three component disk model, the radial profile 
is still sensitive to the structure of the underlying optically 
thick small grain disk surface. In particular, the addition of 
the optically thin component that gives rise to the r -3 power 


law cannot mask a discontinuous small grain distribution at the 
midplane, even when the optically thin component is allowed to 
contribute an order of magnitude more scattered light than the 
optically thick component. A large number of optically thick 
disk models were investigated in this effort, and are described 
in detail in the Appendix. 

5. MODELING DISCUSSION 

Our requirements to fit the radial profile are threefold. First, 
we require that there be enough small grain material inside the 
36 AU large grain cavity to scatter direct NIR starlight before it 
reaches the cavity wall, effectively shadowing it. Although this 
gap must be optically thin to sub-mm light in order to reproduce 
the SMA observations of All, we find no viable scenarios in 
which it can be optically thin at shorter wavelengths and remain 
consistent with the observed H - band radial PI profile. 

Secondly, the smoothness of the observed H - band PI profile 
suggests that the surface of this underlying optically thick small 
grain disk must follow a smooth surface geometry. We cannot fit 
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Figure 16. A schematic view of the three component disk model that provides the best fit to our data. The vertical scale has been inflated by a factor of two relative to 
the horizontal. The origins of sub-mm light, H - band scattered light, and the steep 12-20 /zm SED rise are indicated. 

(A color version of this figure is available in the online journal.) 


our data under scenarios where there is a marked effect, either in 
differential depletion or a scale height discontinuity, of the large 
grain dust gap on the small grain distribution. This is perhaps 
the most interesting and puzzling result of our observations 
and analysis — that the small H - band scattering grains in the 
SR21 disk must follow a smooth distribution that is relatively 
unaffected by the goings-on at the disk midplane. This supports 
scenarios under which dust filtration occurs at the large grain 
gap wall, allowing the small grains to be carried into the gap 
along with the disk’s gas. 

Finally, we find that the steepness of the H - band radial profile 
of SR21 is indicative of an optically thin small grain disk 
component. In order to remain consistent with the first two 
criteria, this requires that we invoke a third vertically extended 
optically thin disk component overlying the optically thick 
portion of the disk. 

As noted in Section 2, the AO system performance was 
sub-optimal during the SR21 observations, and uncorrected 
effects such as those of a polarized seeing halo may mimic 
disk emission, particularly close to the star. Thus, it is possible 
that the true SR21 radial profile flattens out in the interior and 
we are mistaking uncorrected PSF effects for disk emission in 
the inner region. However, all pixels in the saturated region were 
excluded from the radial profiles, and the PSF artifacts do not 
extend beyond the innermost bin in our radial profile plots, so 
this is unlikely. 


The computational expense of scattered light models means 
that we did not investigate the scattered light parameter 
space as systematically as the SED parameter space, and a 
more thorough and systematic exploration of the available 
parameter space is warranted to ensure that our solution is 
unique. 

It could be argued, however, that ADI observations with more 
field rotation, higher AO performance (Strehl ratio) H - band 
observations and/or further multi wavelength data that resolve 
the mysterious SR21 “gap” for other varieties of disk grain are 
equally warranted. Observing methods that are sensitive to sharp 
disk features (e.g., ADI) could help verify whether there is any 
local excess in scattered light at the 36 AU dust wall in addition 
to resolving the potential companion. 

It is also worth noting that all of our disk models have 
assumed axisymmetry. Unresolved spiral structure, emission 
from a companion, disk warps, etc. could be present in the 
SR21 disk and may have an effect on the radial profile that has 
not been examined in this analysis. 

Uniformity of grain composition as a function of radius 
is another aspect of the models that is perhaps unphysical. 
Compositional differences, such as higher volatile abundances 
in the gapped region, could allow for a more marked effect of 
density on surface brightness in the inner disk, allowing optically 
thick large-/! large-a models to increase more steeply toward 
the star. 
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6. CONCLUSION 

We have spatially resolved the SR21 transitional disk for the 
first time at NIR wavelengths. Our data show the signature of 
an optically thin small grain disk component, and place strong 
constraints on the underlying distribution of any optically thick 
surface grains. 

Our Subaru data resolve the disk in scattered polarized 
light for stellocentric r ^ Of! I (>12AU). Extended polarized 
emission is present above noise level out to r ~ O'/ 6 (~80 AU) 
and exhibits a centrosymmetric pattern of polarization vectors. 
The azimuthally averaged radial profile of the PI emission 
reveals that the disk flux drops off steeply (r -3 ) and smoothly 
with radius. The disk is also evident in the PSF-subtracted total 
intensity image, and the range of viable scale factors for PSF 
subtraction give us an estimated polarization percentage for the 
disk of 10% < P < 27%. 

We confirm with an independent model that scenarios in 
which large midplane grains are depleted by a factor of ~10 -6 
inside of a ~36AU gap reproduce the SMA observations of 
All; however, small grains cannot be depleted to the same 
degree inside the cavity. 

Our analysis eliminates a variety of scenarios under which 
there is a discontinuity in the small grain disk at the 36 AU large 
grain gap radius, both in terms of physical likelihood and ability 
to fit the observational data. In order to be consistent with the 
lack of an equivalent gap in the Subaru HiCIAO H - band PI 
images and radial intensity profiles, there must be a decoupling 
of large and small grain disk parameters in the SR21 disk. Our 
data rule out scenarios in which the SR21 disk contains an 
optically thin small grain dust gap with a wall at 36 AU. 

Although we can neither confirm nor deny its existence 
definitively, we present several pieces of anecdotal evidence 
in favor of a companion in the SR21 disk. These include 
an intensity excess near the appropriate location of E09’s 
hypothesized companion, a deviation from centrosymmetry in 
the polarization vectors in this region, and an apparent migration 
of the disk major axis from inner to outer disk, suggesting a disk 
warp or potentially unresolved spiral structure. 

We conclude that the best fitting models for the H - band ra- 
dial PI profile, the observed sub-mm emission and the mul- 
tiwavelength SED for the SR21 disk invoke three separate 
components. These axisymmetric, uniform radial grain com- 
position scenarios suggest an optically thin envelope or atmo- 
sphere component lying above an optically thick, geometrically 
smooth small grain disk component. Further constraints on the 
geometry of this very interesting transitional disk may be ob- 
tained by follow-up at higher spatial resolutions and at other 
wavelengths. 
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APPENDIX 

OPTICALLY THICK SCATTERED LIGHT MODELS 

In this appendix, we detail our attempts to derive a disk model 
in which the H - band scattered light emission originates from the 
surface of an optically thick small grain disk. Although we ulti- 
mately invoked an optically thin component in order to explain 
our data, this was not the most obvious solution, and more 
than 100 optically thick small grain scenarios were investigated 
before we arrived at the conclusion that the steepness of the 
observer H - band radial intensity profile is due to the effect of 
an optically thin small grain disk component. 

Due to the relative expense of full scattered light models, we 
have not explored the entirety of the parameter space for the 
optically thick small grain disk solution. We therefore describe 
our efforts in this appendix, in the hope that they will help inform 
future H - band scattered light modeling efforts. 

A.l. Gapped and Differentially Depleted Models 

If the r <36 AU gap is optically thin at all wavelengths, as 
has been assumed by most previous models, then sufficient flux 
is incident on the dust wall at 36AU to create a sharp rise in 
the SED from 12 to 20 /zm. We call the class of models with 
a depleted cavity that is optically thin in both large and small 
grains “Gapped” models. 

Our initial working assumption, though perhaps naive, was 
that the observed SR21 cavity would also be seen in small grains. 
As discussed in Section 4.3, this hypothesis provides an excel- 
lent fit to the multi wavelength SED. However, if the cavity size 
and depletion factor suggested by the sub-mm data is carried 
over to the small grains that scatter at H- band, the simulated 
H - band PI radial profiles produced by the Whitney code 
(Figure 17) suggest that we should have seen an observable 
discontinuity in the H - band scattered light PI profile at the r = 
36 AU sub-mm wall. We see neither a slope nor a flux discon- 
tinuity at this radius, and can immediately exclude the simple 
scenario in which the large and small grains follow precisely 
the same distribution. 

The simplest scenario that allows us to maintain a similar 
distribution for large and small grains is one in which the two 
populations are both depleted inside of the cavity, but each to a 
different degree (a different 8 value). Physically, this would seem 
to make sense, as the large grains concentrated at the midplane 
should be subject to different physical processes (perhaps most 
notably clearing by forming protoplanets) than the grains at 
the surface of the disk where H - band scattered light originates. 
Second generation small grains are also expected to be generated 
by the planet formation process (e.g., Zuckerman 2001), and so 
even if small grain depletion also occurs at the midplane, it could 
in some sense be compensated for. 

We began by exploring a variety of scenarios in which the 
small grains follow the same basic distribution as the large 
grains, are undepleted in the outer disk, but are depleted to a 
different degree than the large grains inside of the r < 36 AU 
cavity. Radial profiles for models following the same basic 
distribution as the large grain disk (a = 2.15,y0 = 1.15), but with 
depletion values ranging from 10 -8 to zero (10°) depletion, are 
shown overplotted on the observed H - band radial profile and on 
the observed SED in Figure 13. 

In moderate gap depletion cases ( 8 = 10 -3 -10 -5 ), the flux 
drops sharply inside of the gap and falls all the way into the 
saturation radius, which is inconsistent with our observation 
of a continually inwardly increasing profile. In cases where 
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Sample Gapped, Smooth Model 


Large grain density distribution 


Spectral Energy Distribution 



radius (AU) 



Small grain density distribution 


Radial Polarized Intensity Profile 



0.0 33.3 66.7 100.0 

radius (AU) 
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Figure 17. Sample model outputs for a representative gapped disk model. The lefthand panels represent the density of large grains (upper) and small grains (lower). 
The large grain model is identical to the All best-fit model described in Section 4.2. The small grain disk, like the large grain disk, is depleted by a factor of 10 -6 
interior to 36 AU, and it truncated at 7 AU. The righthand panels provide model fits to the observed SED (upper) and H - band radial polarized intensity profile (lower). 
Although the small grain surface layer is continuous under this model, the heavy small grain depletion interior to 36 AU allows a significant amount of flux to be 
incident on the wall, creating a discontinuity in the H - band polarized intensity profile. 

(A color version of this figure is available in the online journal.) 


the depletion is high enough to make the disk optically thin 
(5 = 10 -6 -lCT 8 ), the radial profile well inside the gap does 
increase sharply inward, however there is invariably a brightness 
discontinuity at the gap radius. This is due to the small amount of 
“wall” that is tipped into the line of sight by the disk inclination, 
which creates an excess at this radius that is evident in the radial 
profile. In low gap depletion cases ( 8 = 10°-10 -2 ), the radial 
profiles are relatively smooth, but too flat. 

The “Gapped” models provide excellent fits to the SED and 
sub-mm imagery. However, the increased emitting area of the 
gap “wall” invariably produces a local excess at (f/3. In order 
for such models to fit the observed radial profile, two things 
would need to occur. First, scattered light emission from the 
wall (whose existence is required by the SED fits) needs to be 
somehow suppressed to avoid a sharp discontinuity in scattered 
brightness. If the disk were perfectly face-on, it may be possible 
to hide the excess emitting area of the wall, but the range of 
inclination estimates from a variety of multiwavelength work 
(All; this work; P08) suggest that this is unlikely. 

Second, they require that the optically thin grains inside of 
the gap and optically thick grains outside of the gap share a 


nearly identical radial PI profile despite very different physical 
conditions. Steep radial profiles can be produced by both 
optically thin and optically thick disks, but the apparent precise 
equality between the two makes this scenario seem somewhat 
contrived. We therefore exclude the “Gapped Disk” models as 
incompatible with the scattered light imagery except under the 
most contrived scenarios. 

A. 2. Scale Height Discontinuity Models 

We also consider scenarios that do not require the small grains 
to be depleted (and therefore optically thin) inside of the large 
grain gap. A step-like discontinuity in the small grain scale 
height could contribute to the steep SED rise, without requiring 
that the inner disk be depleted in small grains relative to the 
outer disk and allowing for identical (though discontinuous) 
distributions from inner to outer disk. We call this class of 
models “Scale Height Discontinuity” models, and the outputs 
for a representative model are shown in Figure 18. A related 
class of models utilize a localized “puff” in the large grain outer 
disk wall. Inflated disk rims are often invoked to explain variable 
emission in disks (e.g., Flaherty et al. 2012). Under such models, 
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Sample Scale Height Discontinuity Model 


Large grain density distribution 



radius (AU) 


Spectral Energy Distribution 



Small grain density distribution Radial Polarized Intensity Profile 
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Figure 18. Sample model outputs for a representative undepleted small grain disk with a scale height discontinuity at the 36 AU large grain dust wall. The lefthand 
panels represent the density of large grains (upper) and small grains (lower). The large grain model is identical to the Al 1 best-fit model described in Section 4.2. The 
small grain disk undepleted, with a scale height discontinuity at 36 AU and an inner truncation at 7 AU. The righthand panels provide model fits to the observed SED 
(upper) and H - band radial polarized intensity profile (lower). The large amount of small grain material in the inner disk overproduces the 10 /im flux, but provides an 
otherwise good fit to the SED. In this case the discontinuity in the radial profile at 36 AU is smaller than in the gapped models due to the relatively smaller amount of 
direct stellar flux incident on the small grain dust wall. However, the steep r -3 profile is not reproduced in this case, nor was it with any of the other radial density (a) 
and scale height (/?) exponents investigated for such scale height discontinuity models. 

(A color version of this figure is available in the online journal.) 


the high incident flux on a disk wall causes higher temperatures 
at that location and “puffs up” the rim. If this were to occur at 
the 36 AU large grain disk wall, it is conceivable that it would 
also inflate the small grain scale height at this radius. 

At the limit of a perfectly vertical step-like discontinuity and 
a face-on disk, there is no break in the radial profile between 
inner and outer disk at the H band because the optical depth 
of these grains is large and the wall is obscured by the disk 
surface as viewed from face-on. However, as with the gapped 
disk models, it is difficult to mimic this situation in an inclined 
disk, even one with low inclination like SR21. These models 
thus solve the problem of requiring a different disk geometry 
inside and outside the gap however they also invariably produce 
a local excess at the gap radius that is not seen in the data, just 
as the “Gapped” models did. 

In the case of puffed disk rims, the discontinuity that causes 
the puff is localized to the “wall,” so the disk soon settles back 
into its normal scale height distribution. These “puffs” may 
have a small radial extent, but again it is difficult to construct a 
situation in which they wouldn’t be visible as an excess in the 


radial profile. On the other hand, the puffed disk rim shadows 
a portion of the disk immediately exterior to it from incoming 
stellar flux, where the area of the shadowed region depends on 
the height of the “puff.” This local excess (“puff”) followed 
by a local minimum (shadowed region) could be smoothed 
out in the PI radial profile if the total radial extent is much 
smaller than the diffraction limit at H- band, mimicking our 
observed smooth radial profile. This precise balance again seems 
somewhat contrived, so we consider the puffed rim explanation 
to be relatively unlikely, and we do not include any of these 
models in our final fits. 

Owing to their continuous density distribution, brightness 
discontinuities in the “Scale Height Discontinuity” models are 
less drastic than those created by the “Gapped” disk models, so 
they are more likely to remain unresolved in our observational 
data. Non-diffraction limited performance, the peculiar shape 
of the HiCIAO PSF, and the relative insensitivity of PI emission 
to sharp disk features (due mainly to low overall flux) could all 
conceivably smooth out such discontinuities in the radial profile 
to a greater extent than convolution with a (T06 Gaussian PSF 
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Sample Steeply Curved Model 


Large grain density distribution Spectral Energy Distribution 
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Small grain density distribution 
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Figure 19. Sample model outputs for a representative model in which the region of the sub-mm cavity is more steeply curved in the small grain disk. The lefthand 
panels represent the density of large grains (upper) and small grains (lower). The large grain model is identical to the All best-fit model described in Section 4.2. 
The small grain model is undepleted interior to 36 AU and truncated at 7 AU. The scale height exponent /3 was increased from 1.15 in the outer disk to 2.5 interior to 
36 AU, allowing for a steeper curvature. The righthand panels provide model fits to the observed SED (upper) and H - band radial polarized intensity profile (lower). 
This small grain distribution completely masks any discontinuity at the 36 AU dust wall created by an H - band contribution from the large grain dust prescription; 
however, the steep inner disk carries a flat radial profile, which is inconsistent with the observed r -3 profile. 

(A color version of this figure is available in the online journal.) 


would suggest. However, we were unable to construct a disk 
model in which the radial profile was smooth enough to be 
consistent with our data and error bars, so we exclude these 
models as relatively unlikely candidates. 

A3. Steeply Curved Models 

If we relax the assumption that the small grain disk needs to 
have precisely the same parameters inside and outside the gap, 
we find that there are cases where a turnover in the radial scale 
height exponent /3 at the 36AU wall can reproduce the steep 
SED rise. This situation does not require small grain depletion, 
nor does it require a discontinuity at the large grain dust wall. In 
such cases, the SED rise coincides with the region just inside of 
36 AU where the disk is most steeply curved, and in that sense 
mimics a “wall.” From a physical standpoint, this transition from 
a large to a smaller radial scale height exponent at the large grain 
dust wall would need to be explained somehow by the dearth of 
large grains at the disk midplane in this region or the physical 
mechanism causing disk clearing, which we do not endeavor 
to do here. We call this class of models “Steeply Curved” 


models, and outputs for a representative model are shown in 
Figure 19. 

Under the third well-fitting SED scenario, a steep disk 
curvature (high /3) in the inner disk mimics a “wall” at the 
transition to a shallower outer disk profile. Without a physical 
discontinuity in the small grain surface at 36 AU, these models 
are also free of discontinuity in their radial PI profiles. However, 
all “Steeply Curved” models also produce flat (r°) radial profiles, 
which our observations argue against. This is because the 
geometric emitting area per pixel increases as the disk becomes 
increasingly flared, so pixels represent a larger amount of disk 
surface as you move outward. This mitigates the r -2 effect 
of incident stellar radiation, and makes the flux in the disk 
extremely flat. 

We find that the surface brightness of the disk is relatively 
insensitive to the radial density exponent a in our high-/3 models, 
even for very strongly inwardly peaked density profiles (e.g., 
r -10 ). These models still show a flattening of the brightness 
profile in steeply curved regions that our data rule out. We 
therefore exclude these high-/3 undepleted scenarios as poor fits 
to the H - band data. 
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Table 3 

SIRPOL Aperture Polarimetry 


Band 

Polarization 

PA 


(%) 

(°) 

J 

2.32 ±0.17 

21.5 ±2.1 

H 

1.01 ±0.05 

18.0 ± 1.3 

Ks 

0.52 ±0.10 

179.2 ±5.3 


We find that the steepness of the radial profile significantly 
limits the disk geometry, particularly in the optically thick case. 
We tried more than 30 iterations of smooth optically thick disks 
with a wide range of a, f and a-f> (the radial dependence of the 
surface density) values, and were unable to create a steep r -3 
profile. 

The only optically thick disks we found that reproduce an r~ 3 
radial PI profile are geometrically flat. As the geometrically flat 
case requires that the scale height of the disk be smaller than 
the radius of the star, it is impossible to create one in which the 
large grain disk required to reproduce the sub-mm flux does not 
stick up above the small grain disk. It is difficult to imagine a 
situation that would reproduce this, so we exclude these models 
on physical grounds. 

A. 4. Prospects for Envelope Component 

Since SR21 may still be weakly accreting (the upper limit 
of 10 -8,84 on the accretion rate is relatively high compared to 
other disks), we also explored the possibility that an envelope 
component is contributing to the flux in the interior of the gap. To 
provide observational constraints on the amount of circumstellar 
material, SR21 was observed on 2011 March 29 with SIRPOL 
(Kandori et al. 2006), a ///^-simultaneous imaging polarimeter 
that is mounted on the Infrared Survey Facility (IRSF) 1.4 m 
telescope in South Africa. We conducted aperture polarimetry 
on this data and the results are summarized in Table 3. Our data 
reveal total polarization levels that are consistent with the range 
of total polarization observed for YSOs (Pereyra et al. 2009). 
We have not corrected the SR21 aperture polarization values 
reported in Table 3 for interstellar polarization, which could be 
significant due to the high foreground extinction ( A v = 6.3). 
Moreover, because the observations were conducted in defocus 
mode to avoid saturation, the A and B components of SR21 
(separated by 6") overlap in our aperture polarimetry, preventing 
us from using SR21B to determine how much of the measured 
polarization is due to foreground dust and how much is local to 
SR21A and/or B. There are no other nearby objects in the field 
suitable for use as field star interstellar polarization probes. 

If we assume that the level of interstellar polarization and 
contamination from SR21B is minimal, we can use the measured 
polarization Position Angles (PAs) at each band to probe of 
the optical depth of the disk, as optically thick disks show 
polarization PAs oriented perpendicular to the PA of the disk 
major axis, while the polarization PA and major axis PAs of 
optically thin disks are oriented parallel to one another (Whitney 
& Hartmann 1992; Pereyra et al. 2009). The disk major axis PA 
(86° ± 1 1°, this work; 100°, Al 1) is roughly perpendicular to the 
measured polarization orientations at all three bands, suggesting 
that the disk may be optically thick in the NIR. 

However, there are several caveats to be made. First, the 42° 
counterclockwise migration in the disk PA derived from isopho- 
tal fits (and described in Section 3.1) suggest that an average 
PA is not a particularly good measure of the overall disk orien- 
tation. Secondly, these results depend on the robustness of our 


assumption that the interstellar polarization and contamination 
from SR21B is small. Finally, and most importantly, we have 
managed to eliminate all of the optically thick small grain dust 
disk scenarios that we have investigated thus far. In fact, we find 
that optically thin disks provide a much better fit to the H- band 
radial profile, as described in Section 4.4, although we cannot 
rule out, and in fact require, an additional contribution from an 
underlying optically thick small grain disk component. 
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